Localised projective measurement of a quantum field in non-inertial frames 
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We propose a projective operator formalism that is well-suited to study the correlations of quan- 
tum fields in non-inertial frames. We generalise a Glauber model of detection of a single localised 
field mode that is capable of making measurements in an arbitrary reference frame. We show that 
the model correctly reproduces the Unruh temperature formula of a single accelerated detector, and 
use it to extract vacuum entanglement by a pair of counter-accelerating detectors. This latter ex- 
ample is a proof of principle that this approach will be appropriate to further studies on the nature 
of entanglement in non-inertial frames and, in general, to model experimentally feasible scenarios 
in quantum field theory in non-inertial frames. Finally, as further confirmation of the validity of 
our approach, we introduce an explicit perturbative matter-radiation interaction model which re- 
produces both the generalised Glauber model and the projective measurement results in the weak 
coupling regime. 
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I. INTRODUCTION 

One of the most fascinating features of relativistic 
quantum fields is that the notion of a particle is frame 
dependent. Even if a given observer describes the field 
state as free from particles, the state of the same field de- 
scribed by another observer can be in fact populated by 
particles. This relativistic effect is a direct consequence 
of the Bogolyubov transformation between inertial and 
non-inertial reference frames and has inspired the devel- 
opment of the whole new field of relativistic quantum in- 
formation. The main objective of this study is to find the 
consequences of relativistic effects on protocols describing 
the storage, processing and transmission of quantum in- 
formation. Perhaps most striking of these investigations 
has been the realisation that entanglement, which is a 
key resource in quantum information protocols, depends 
on the state of motion of the observer. 

Seminal works on this topic [1 studied the effects of 
acceleration on the entanglement of global states of free 
fields. However, several issues remained unsettled, in- 
cluding how to properly describe a possible experimental 
setting [2 and how to correctly address the technical so- 
lutions beyond the single mode approximation |3]. The 
recent approach ^ reanalysed the general setting of [1] 
beyond the single mode approximation. However, it re- 
quired the preparation of a whole family of orthogonal 
states of global Unruh modes, one for each accelerated 
observer, instead of a single, fixed state given to all ob- 
servers. It is was also not clear how such global modes 
could be physically prepared. For these reasons the phys- 
ical effects of acceleration on entanglement remained in 
disrepute. 

In this paper we describe a framework for which these 
issues can be resolved. We introduce a model of a detec- 
tor that performs a projective measurement on a localised 



single- mode of the quantum field. Our detector can ei- 
ther move inertially or accelerate. This framework allows 
a fixed state of the field defined in an inertial frame to 
be probed for any acceleration. Thereby avoiding having 
to use a different initial state for each acceleration and a 
non-local initial mode. 

To motivate the projective measurement formalism 
we introduce a relativistic generalisation of the Glauber 
model U] and show the relationship between the quanti- 
ties calculated in each approach. One particularly nice 
feature of the Glauber detector is that the detector only 
clicks when particles are present in the mode that the 
device measures. This lends itself to a particle detector 
interpretation rather than a detector of fluctuations like 
that of the Unruh-DeWitt detector [5] . 

The model is applied to the simplest case of the 
Minkowski vacuum state and we use it to reproduce the 
well-known Unruh effect for a single accelerated detec- 
tor [6J . We then study how a pair of such detectors can 
extract nonlocal correlations present in the vacuum [7], 
showing that our model is suitable for studies of entan- 
glement in non-inertial frames. Our analysis can be ex- 
tended to arbitrary states of the field, which is the subject 
of another work [5] . 

We can freely choose the size and the shape of the 
detected wave-packet mode allowing one to study how 
these properties affect the measurement results. We find 
that the entanglement extracted from the vacuum by a 
pair of uniformly counter-accelerating detectors increases 
with the size of the wave-packet mode and the detector's 
acceleration. 

We first outline, in section |IT] the conceptual ideas be- 
hind the detector model presented here. Then, in section 



III we introduce the formalism to work with uniformly 



accelerated observers of a quantum field and discuss the 
process of field quantisation in the Rindler basis. We give 




FIG. 1. Detection of a quantum state of a single light-mode 
by means of a lens L, single mode fibre SMF and a photon- 
counting detector D. 



then a heuristic argument for a natural choice of physi- 
cally meaningful localised detector modes in section |IV[ 
In section W\ we calculate the average number of parti- 
cles populating the field mode probed by a localised ac- 
celerated detector when the field state is the Minkowski 
vacuum state; and in section [VT] we introduce the formal- 
ism for the projective measurements on the field mode 
probed by the detector. This will allow us to estimate 
the temperature of the photo-count statistics (in section 
VII ) and to study the extraction of field entanglement of 



spatially separated spacetime regions (in section VIII I 



Our conclusions are presented in section |IX| Addition- 
ally, a detailed analysis on the photo-detector interaction 
dynamics is given in appendix \K\ 



II. MODEL OF THE DETECTOR 

We consider a model of measurement carried out lo- 
cally on a quantum field with a device that is only sensi- 
tive to a single localised wave-packet mode. The model is 
inspired by a quantum-optical setting depicted in Fig. [T] 
that has been considered both in theory [9^ and realised 
experimentally [10^. The original scheme consists of a 
lens, L, focusing incident light onto one end of a single- 
mode optical fibre SMF that transmits only one par- 
ticular transverse mode of light. The other end of the 
fibre is attached to a single-photon detector, D, utilising 
the photoelectric effect to detect and count single light 
quanta arriving in the fibre-coupled mode of light. In 
the idealised version of this scheme the state of the sin- 
gle field-mode is projected onto the Fock basis and all 
the other modes that are orthogonal to the one transmit- 
ted through the fibre are not affected by the detection 
process. This scheme can also be generalised to allow 
the detection of an arbitrary mode by placing optical el- 
ements in front of the lens, L, which transform the given 
input mode into the mode transmitted by the fibre. 

Notice that the state of the measured system is com- 
pletely destroyed during the measurement process (pho- 
tons are absorbed) and the physical measurement corre- 
sponding to this model is not instantaneous. The mini- 
mum time of such a detection is equal to the time it takes 
light collected by the lens to reach the SMF plus the time 
it takes for the device that captures the light to make 
the read-out. Therefore the bigger the mode, the longer 



it takes to extract it from all the other modes. In many 
practical scenarios, such as the one considered above, the 
time of detection can be safely neglected. However, one 
has to bear in mind that such idealisation has its limits 
and not taking measurement time into account can lead 
to problems |11| . 

In this work we consider an abstract "black-box" de- 
vice that measures a wave-packet mode of the field in 
the same way: the process consists of filtering out the 
relevant mode and then measuring its state. For simplic- 
ity we consider a 1+1 dimensional version of the above 
setting. The process we describe disturbs the field only 
in the proximity of the measurement device; we assume 
that the detector's response is not affected by the state 
of other modes orthogonal to the one it measures, and 
that the measurements performed on this mode do not 
affect the state of the orthogonal modes. In this sense 
the process is localised in space. 

The quantum field under consideration is taken to be 
a real non-interacting massless scalar field, 0, described 
by a Klein-Gordon equation, Dc^ = 0, in 1+1 Minkowski 
spacetime. Let the wave-packet observed by the detector 
be denoted by lij-oiS., t), where ^ and t, are coordinates of 
the reference frame co-moving with the detector. We as- 
sume that ^d(Si ''") is consists only of positive frequency 
waves with respect to r. Although the splitting into pos- 
itive and negative frequencies is frame dependent, the 
linear combination ipr>{£,,T) + iP^{£,,t) does not depend 
on the coordinate system. It is this real field solution 
that describes the classical wave which couples to the 
detector. Using the positive frequency part it is possi- 
ble to define an operator d = (f/'D,'^), associated with 
the wave-packet that the detector sees. Here (•, •) is the 
Klein-Gordon scalar product defined on the set of com- 
plex solutions to the Klein-Gordon equation by: 



(01,^2) =« 
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g(j)l{x,t)dtMx,t)\t=o, (1) 



tr> 



where g is the determinant of the metric and dt is 
the two-edged differential operator defined by: f dth — 
f c5i~ cdi^' ^^ '^D is unit normalised with respect to the 
Klein-Gordon scalar product, then d satisfies the usual 
commutation relations [d, d^] = 1. If the motion of the 
detector is inertial (uniformly accelerated) then d annihi- 
lates the Minkowski (Rindler) vacuum state. This leads 
to the identification of d'^^^ as an annihilation (creation) 
operator. Any state of the field that is formed by the 
linear superposition of the basis states {-^^=^10)} or a mix- 
ture of them is said to be a state in the mode tpD . The 
operator d' d is a number operator for the number of par- 
ticles in the mode V'D- These excitations are, however, 
not energy eigenstates of the free-field hamiltonian. The 
state accessible to the detector is the total field state |^) 
projected onto the ip-Q subspace: 



Trifel*)(*l 



(2) 



where the trace is to be taken over the subspace of all 
modes orthogonal to '0D- In this work we assume that 
tpD is localised [T5], which will allow us to attribute the 
detector a certain position in space. For example if the 
measurement on ipD takes place at times close to r = 
then the detector has to be placed in the region where 
■0d(^,O) is non-zero. 

One of the greatest advancements of the field of 
quantum optics occurred with the introduction of the 
Glauber model [3]. The detection process described in 
the Minkowski frame by the Glauber model for the scalar 
field results in a particle absorption occurring at the point 
{ct,x), described in the weak coupling regime by the fol- 
lowing absorption amplitude: 



Afix,t)^{f\^+ix,t)\i), 



(3) 



where \i) and |/) are initial and final states of the field 
and (j)^{x, t) is the positive frequency part of the field op- 
erator in Minkowski coordinates. However, the final state 
of the field remains unknown after the measurement, so 
the total probability V of observing a particle is given by 
taking the squared modulus of the amplitude and sum- 
ming over all final states. Generalising to mixed states 
and substituting ^ for the initial state results in: 

P(x,t)=Tr{<^-(x,t)0+(a;,t)TriV*l*)(*l}- (4) 



In this work we generalise the Glauber theory of de- 
tection to the case of detectors moving with relativis- 
tic accelerations. The crucial step in this generalisation 
is to replace the positive and negative frequency parts 
in the above equations with the splitting with respect 
to the proper time of the frame co-moving with the de- 
tector {cT,^). Since the positive and negative frequency 
parts of the field are different in the accelerating and in- 
ertial frames, the intensity of particle clicks in the photo- 
counter will depend on the detector's motion. We prove 
in appendix [X] that our generalised Glauber model agrees 
to first order in perturbation theory with the predic- 
tions of a dynamical field-apparatus model that utilises 
an Unruh-DeWitt interaction, analogously to the inertial 
case [T2]. We find that even in the non-inertial case the 
measurement outcomes can be completely described by a 
set of projective operators, as in the standard quantum- 
optical non-relativistic setups. 




FIG. 2. Minkowski diagram, with Rindler regions I (s > c\t\) 
and II (x < —c\t\) covered with conformal Rindler coordinates 
(r, ^) and (t',^') respectively. 



erated observers are the Rindler coordinates: 



rf=-e</=%inh- 
a c 



ct=^e<'/^%inh — 
a c 



X = —e""^!"' cosh — , 



(/) 



ar 



-e°«'/^ cosh— , (//), 
a c 

(5) 



where — oo < ^, ^', t, t' < oo, a is an arbitrary parameter, 
and /, // represent the two distinct regions of spacetime 
shown in Fig. [2] To interpret a as the proper acceleration 
of the detector, we must ensure that the photo-counter 
is spatially localised around ^ = when a > 0, or ^' = 
when a < 0. Then all of its components will approxima- 
tively experience the same proper time, which coincides 
with Rindler time t. 

One can decompose the field operator into a countable 
sum of Rindler modes by putting the field in each re- 
gion into a box of length h satisfying periodic boundary 
conditions. Then: 

= ^ bk,iWk,i{(., t) + h.iiWk.iiiS,', t') + h.c, (6) 

k 



III. THE ACCELERATED FRAME OF 
REFERENCE 



In this paper we focus exclusively on measurements 
made by the uniformly and relativistically accelerated 
detectors. A natural choice of coordinates for the accel- 



where fc = ^ with n = ±1, ±2 

h ' 






and 

;*(fcC-|fc|cr)g(^_^|^|)^ (7) 
,^(-fe?'-|fe|cr')g,(_^. _^|^|)_ 



The Rindler modes are box normalised so that 
(wkj,Wk'j) = Skk'- We will also find it useful to re- 
express our results in the continuum limit h ^ oo. This 
is done in the usual way by replacing ^ X]fc ~^ I ^^ 

and \ h^Wkj —^ Wkj where in the continuum limit the 
Rindler modes are plane-wave normalised (wkjjWk'j) — 
5{k — k'). The region I and II Rindler modes are posi- 
tive frequency with respect to the Killing vectors idr and 
idr' in region I and II respectively, and have been chosen 
such that fc > corresponds to a right moving wave; the 
asymmetry in sign between equations ([7| occurs because 
^ increases in the right direction while ^increases toward 
the left. 

It is a known fact |i^ that the Minkowski vacuum 
state can be expressed as a multi-frequency two-mode- 
squeezing of the Rindler vacuum state: 



|0)m = Siji\0)r, 



(8) 



where the squeezing operator, Sjji, is characterised by 
the squeezing parameter, r^ — arctanh(e~'^l*''l^ /°), and 
fulfils the following relations: 



S\jjbijSiji = coshnS;^ + sinhr;6j_^^. 



(9) 



We will find these commutation relations extremely use- 
ful when performing the algebraic manipulations of the 
following sections. 



IV. CHOICE OF THE DETECTOR MODE 

The wave-packet corresponding to the localised field 
mode probed by the detector iPd{S,,t) must be spatially 
localised around ^ = at the measurement time around 
r = 0. In principle, any positive frequency localised 
wave-packet satisfying this condition will be sufficient, 
but our choice will be dictated by the convenience of fur- 
ther calculations. 

Consider a resting detector that selects ip-oix^t) such 
that it has a central frequency coinciding with the A^- 
th eigenmode of a cavity. We ask: what wave-packet 
tjJD{^,T) would a uniformly accelerated detector select 
if it was to operate on this same principle? In order 
to answer this question we compare the eigenmodes of 
an ideal cavity at rest and the eigenmodes of the same 
cavity in uniform acceleration [13!. One infers that the 
central frequency of the mode tlju{^,T) should be — , 

where, a = ^asinh (^), L is the proper length of the 
cavity and A^ is a natural number corresponding to the 
chosen cavity eigenmode. We choose a Gaussian envelope 
of that mode, characterised by the width ct, to make it 
localised, however all our main results will not depend 
on the specific choice of that envelope: 



Mto) 



drMi,0) 



1 



^/n\/2^ 



cxp 



-T^+*-f 



Nc 

-i — 7/.D(e,0), 



The above expressions always have some contribution 
from negative frequency components, however these are 
negligible for A^ ^ 1 (actually A^ > 3 is already enough) . 
To get rid of them completely we introduce an infra-red 
cut-off A to the considered mode: 



V'D(e,^) 






(11) 



Note this then also ensures that the single particle states 
exist for this degree of localisation [15 . The device is 
defined to be only sensitive to the right-moving modes in 
region I but this choice is arbitrary and can be changed 
if necessary. In the limit of small accelerations, the 



mode ( 10 1 written in Minkowski coordinates reduces to 

1 



V'd(2:,0) 



exp 



-^2 +i L^ 



up to a translation 



along X. It reproduces a regularised A^-th mode of a rest- 
ing cavity of a length L. 



V. DETECTION OF PARTICLES FROM THE 
VACUUM STATE 

Consider measuring the occupation number operator 
of the transmitted localised mode, (Pd. If at r = 
the measured field is in a Fock state \n) of a mode 
(J){£,,t) that coincides at r = with V'd, i-e., 0(^,0) = 
^d(^, 0), 9t</'(^,0) = 9r'0D(C;O) then the outcome of 
the detection will be a classical variable n. 

To generalise the result to an arbitrary state of the 
field it is useful to introduce the following decomposition 
of the field operator: 

0(e, 0) = diPuiL 0) + dVolC, 0) + 0', (12) 

where 0' is the remaining part of the mode decomposi- 
tion containing all the modes orthogonal to ipD. In prin- 
ciple, a complete and orthonormal basis could be found 
for the Hilbert space by using Gram-Schmidt orthonor- 
malisation on the positive frequency subspace of complex 
solutions beginning with the vector ^q. However, it will 
not be required to calculate such a basis for anything that 
follows. 

The probability V of the detector clicking at time r = 
and position ^ according to the generalised Glauber 
model Q is: 



^(e,0) = |VD(C,0)|2(M'|dtd>) 



(13) 



(10) 



Of course, the Glauber model describes an idealised sit- 
uation in which the measurement is performed instan- 
taneously and at a localised position. It also ignores 
the fact that the apparatus itself is a complex (classi- 
cal and quantum) system. Some of these considerations 
are taken into account in Appendix |^ where a model of 
field-apparatus interaction is considered. It is neverthe- 
less found that up to an apparatus dependent propor- 
tionality constant (i.e., dependent on collection efficien- 
cies etc) the single click probability is still proportional 



to the expectation of the operator (vd, i.e., the number 
of particles in the mode. 

Let us calculate the expectation value of the number of 
particles seen by an accelerating detector when the field is 
initialised into the Minkowski vacuum state, j'l') = |0)m: 



Sd 



iru,m^iy,$) 



(14) 



From here on the angular brackets represent taking ex- 
pectation values with respect to the Minkowski vacuum 
state. Using (|8| and (l9| one finds: 

k,l 



X K{0\SljjbljkjSiji\0)K 

= ^{nk)\{ipD,'Wkj)\'^, 



(15) 



where {fik) = sinh r^ is the average occupation number 

of a two-mode squeezed state of the plane-wave mode k. 

Let us calculate the number of particles for the mode 

(10 1. Since we will perform several calculations in this 



paper using the same technique, we will go through this 

2 

calculation in detail. We consider the regime, fj <C 1, 
so (T ~ L. If we take the cut-off at A ^ 1/L and make 
the change of variable k' — kL we find: 



(d'td) 



■-{N- 



ANJ2ti 



dfc'(7V + fc')'e-5(fc'-^-^)' 



1 



(16) 



where we completed the square in the exponent, and took 
out the fc'-independent factor. The denominator of the 
second factor in the integrand is approximately unity on 

9 2 

the whole domain, furthermore in the limit N ^ ^^^ the 
integral is dominated by the contribution from k' ^ N. 
Thus, we find up to an order 1 proportionality constant: 



(^d^d) 



-(N^ 



'-) 



(17) 



This equation is only valid for the range of parameters 
satisfying the conditions above, namely 5^ 3> ^ ^ 1. 
As we expect, the number of particles that can be ob- 
served by the detector is exponentially suppressed by the 
factor e~2'^'^"=/'' where Wc = ^ is the characteristic fre- 
quency of the detected mode. 



VI. PROJECTIVE MEASUREMENT 
FORMALISM 

Notice that calculating the expectation value of the 
number of particles in the mode can be done in a purely 
abstract way without any reference to the detection pro- 
cess. In fact, working with the field theory at this ab- 
stract level of description can often be simpler and more 



illuminating than getting lost in the details of a specific 
detector model. The aim of this paper is to show that 
the usual phenomena of fields in non-inertial frames can 
be re-derived at this abstract level. 

The notion of a quantum state makes sense, because 
in principle there could exist a device that measures it 
acquiring the complete information about the system en- 
coded in the state. Appendix [^ shows that the abstract 
calculations of the previous section can also be obtained 
with a dynamical model of interaction between a small 
detector and a quantum field. In a similar way all the 
information extract able from a device that resolves the 
number of particles in a state must be describable by a 
set of projective operators \n){n\ associated with possi- 
ble measured outcomes. This is imposed by the Hilbert 
space structure of quantum states. 

To calculate the probability 7^(n) of detecting n field 
quanta we use the projector onto the n-th Fock state of 
the mode d , that can be written in the following manner 



V{n) 



-d'*d 



(Sdy 



(18) 



In the weakly coupled regime the above expression can 
be derived from a dynamical model described in the Ap- 
pendix LAI with an extra modification of dW being re- 
placed with rjd^d, where rj < 1 |12| . The new parameter 
T] describes the efficiency of the detection and is related 
to the coupling strength between the systems. 

The probability distribution ( 18 1 can be related to the 



characteristic function Z{X) via the Fourier transform: 



dA 

— f 

27r 



-iXn 



Z{\). 



V{n) = 

Here 

00 
Z(A) = ^e*^"7'(n) 



(19) 



= (e'^-l)dt 



E 



(e' 



^ sinhrfc(V'D, Wkj)*bkJi 



y^sinhn(VjD,w;j)6| 



// 



L ; 



= ^(e^^-l)"(dtd)n^_ 



(e»^-l)(rftdy 



(20) 



where the multinomial expansion was used to obtain the 
last line and the mean particle number (d^d) is given by 
Eq. (15). Through (19 1 we find the excitation statistics: 



{Sdy 



V{n) = 

(l + (dtd))l+" 

for the Minkowski vacuum state |0)m- 



(21) 



VII. DETERMINING THE UNRUH 
TEMPERATURE 

One might wonder what temperature is associated with 
the above statistics. For a thermalised state of a har- 
monic oscillator the temperature is defined by the rela- 
tion (d^d) — [e^^^^ — l) , where E is the energy of 
a single excitation of the field mode. However in gen- 
eral the mode i/jd has a frequency spread so it does not 
have well defined energy, therefore the temperature is 
also not well defined. Nevertheless, one can define a tem- 
perature estimator Tost by replacing E with the expec- 
tation value of the energy of a single excitation of the 
field mode, E -^ R(0|(ii^Rdt|0)R = ^^ /k^feKV-o, w'/c,/)!', 

where Hn = J2k f^kbk i^kj + {I ^ II) is the gener- 
ator of proper time translations in the Rindler frame 
and LUk are Rindler frequencies. For a detector moving 
with proper acceleration, a, this temperature estimator 
is equal to: 



r^Jfist — 



X]fcfe/^/c|(V'D,Wfc,/)|' 

log(l + (dtd)-') 



(22) 



Here we find that the temperature observed depends on 
the proper acceleration of the detector, a, and the shape 
of the mode tpD . 



For the choice of mode ( 10 1 , we find that the tem- 



perature estimator ( 22 1 reduces approximately (when 



N 



» h » 1) to: 



/CJ o' 



ha 



2-Kcl~TTC^IaLN' 



(23) 



Note that the temperature estimator deviates from the 
Unruh temperature formula jH] for finite N . This stems 
from the fact that the energy of the mode is not well- 
defined. However, for large N , which corresponds to a 
peaked energy spectrum, we recover the celebrated Un- 
ruh result [6J. 



VIII. NON-LOCALITY OF THE VACUUM 
STATE 

In order to extract entanglement from the vacuum 
state we need more than just one detector. Consider 
two identical detectors, one moving with proper accel- 
eration a > in region / and coupled to a mode ^d,/- 
and the other moving with acceleration —a in region // 
and coupled to ipB.ii- The corresponding annihilation 
operators are dj and dn respectively and both detec- 
tors are causally disconnected. One can generalise the 
calculation ( 15 ) to determine the average product of the 



particle content detected in the two modes and obtain 
the result: 



djdidjjdii 



d\jdii 



djdn 



(24) 
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FIG. 3. Entanglement estimator £ as a function of the di- 
mensionless parameter ^ ior N = 800 (solid line), A*' = 1200 
(dashed line) and iV = 1600 (dotted line). 



where the first two terms on the RHS of the equation are 
given by ( 15 ) and the last term is equal to: 



didii) ^ Yl \/{nk){l + {nk)) 



x{-iprij,Wkj)iipB,ii,Wkji). (25) 

We find that the measurement outcomes carried out by 
the detectors are correlated and the correlations are de- 



termined by the magnitude of expression ( 25 1 . In order 



to prove the non-locality of these correlations we note 
that a two party symmetric Gaussian state is entangled 

ifUZl: 



((i/ - xiif^ {{pi+Pii)j < 1, 



(26) 



/2x„ 



do- + dl and ^/2ipa — da — d\, for a G 



where a 
{/,//}. 

In order to detect entanglement we need to re-program 
our detectors to perform projective measurements cor- 
responding to the hermitian quadrature operators x„ 
or Pa instead of projecting in the Fock basis. Typi- 
cally, such a measurement is realised by means of ho- 
modyne detection with the use of an auxiliary beam of 
light [18'. In the present work we do not study the de- 
tails of such setup [19^ assuming that our device carries 
out an ideal measurement of the quadratures. We find 
that the left-hand side of the inequality ( 26 1 is equal 



to 



( 1 -I- ( d\di ) + ( d\jdji \ — 2Re ( d/d// ) j and again 
it turns out that the presence of the entanglement is dic- 
tated by the same term (25 1 that was responsible for the 



existence of correlations between the detectors' counts. 
For a pair of identical, counter-accelerating detectors we 

have (djdj) = (djjdjj). By this criterion we have 

checked that for any finite acceleration the state mea- 
sured by the detectors is entangled. Knowing this, we 



can now quantify this entanglement extracted by the de- 
tectors by introducing the entanglement estimator £: 



£ = log 



Re ( djdjj 



C 



(27) 



where C is an arbitrary real constant factor. Here 
£ — > — oo when the detectors state is separable and grows 
monotonically as the degree of quantum entanglement is 
increased. For V'd./ given by (10 1 and V'd.// given by the 
transformation ^' = — ^ of (10 1 we have {'ipujiT'u^kji) = 

{ipujjWk.i)- Taking the limit of 
mator can be approximated by: 



2. » II » 1 tlie esti- 



£ 



TTC 

aL 



N - 



2aLj 



C. 



(28) 



We know that the extracted entanglement must vanish 
{£ — > — oo) as a — )■ 0, although for fixed N our approx- 
imations are not valid in this limit. In Fig. |3] we plot 
the dependence of the entanglement £ ( p7| on the di- 
mensionless quantity ^ for several values of N . find- 



ing good agreement with the approximated result ( 28 1 . 
We find that the largest amount of entanglement can be 
extracted from the Minkowski vacuum state when the 
proper lengths and proper accelerations of the detectors 
are large and the frequency numbers N are low. 



IX. DISCUSSION AND CONCLUSIONS 



Although it is both mathematically and conceptually 
advantageous to work at the abstract level of projective 
operators, one must keep in mind the underlying physical 
processes. For example, the state of the measured system 
is completely destroyed during the measurement process 
(photons are absorbed) and the physical measurement 
is not instantaneous and not taking them into account 
can lead to problems |11| . In Appendix E we showed by 
taking the interaction time into account that under rea- 
sonable conditions the results of the generalised Glauber 
model can be validated. 

The model presented here can be readily applied to 
other types of motion, as well as curved spacetimes in 
regions where time-like Killing vectors exist. Further re- 
search is currently being undertaken to use this approach 
to investigate localised entangled states, instead of the 
Minkowski vacuum. The principles laid out in this work 
can be generalised to other types of fields, and it is ex- 
pected that new insights into the nature of entanglement 
in accelerated frames are to be found in this approach. 
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Our results may help to understand the nature of the 
correlations between localised measurements of the quan- 
tum vacuum studied in the context of algebraic QFT 
|20| . It is also interesting to compare our findings with 
the results of pi], where entanglement between two lo- 
calised, orthogonal modes was observed in the inertial 
frame. This entanglement was shown to decay with the 
spatial separation of the modes. On the other hand, we 
find that no entanglement can be extracted from inertial 
detectors. This difference can be attributed to the fact 
that our detectors are probing modes that consist only 
of positive frequencies and the authors of [21_ consider 
modes that have no such restriction. 

The main idea of this paper is to emphasise that for 
single mode measurements the outcomes of photo-count 
experiments in the accelerating frame can be described 
by projective operators operating only on the single mode 
subspace. In this way calculations can be performed at an 
abstract level without making any reference to the actual 
detection processes. The versatility and power of this 
approach is demonstrated by the analytic expressions we 
were able to re-derive for the standard phenomena in non- 
inertial frames i.e., the Unruh effect and extraction of 
vacuum entanglement. Although here we restricted our 
attention to measurements of the occupation number of 
the mode (see Section Q and quadrature measurements 
(see Section VIII ) , in principle any observable in question 



Appendix A 

In this appendix we consider a dynamical interaction 
model for the photon absorption processes occurring at D 
in Fig. [l] and show that the model reproduces the general 
behaviour of a generalised Glauber model (13 1. 



In idealised photo-detectors electrons in bound po- 
tential wells are excited into free states by the absorp- 
tion of photons, see for example [12!. In this approach 
the interaction of the photon with the electron is cal- 
culated using a non-relativistic Schrodinger equation for 
the electron and using the minimal coupling substitution 
p — > p — eA{x, t) resulting in an interaction term. 



Hi = —p- A(rQ,t). 
m 



(Al) 



Since the transformation of this setup into the acceler- 
ated frame is complicated by many factors (including rel- 
ativistic considerations, and the effects of acceleration on 
the binding material properties etc) it is more practical 
to consider a detector with subsystems interacting with 
the field via the interaction Lagrangian: 



L[x{t)\ = p.m{T)(j)[x{T)] 



(A2) 



can be chosen at will by the experimentalist. 



where x{t) is the classical path followed by the detec- 
tor, T is the proper time along the path and to is a 
monopole moment having discrete energy levels. The 
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time evolution of the monopole moment can be writ- 
ten, m{T) = e*^o'"/'*TO(0)e"*^«"'/'* where Hq is the free 
Hamiltonian for the monopole. This is of course just 
an Unruh-DeWitt interaction j21]. For an inertial path 
equations (All and (A2| are remarkably similar. We can 



even approximate the continuum of energy states in the 
photo-electric case (above the binding energy > hcuo) by 
discrete energies of the monopole. For example, we could 
assume that the energies of the monopole have a large gap 
from ground to first excited state representing the bind- 
ing energy of the electron potential then smaller perhaps 
equally spaced gaps thereafter. One would then be able 
to reproduce the standard theory of photo-detection us- 
ing this model |12| . Furthermore, since the interaction 
transforms simply as a scalar under coordinate transfor- 
mations, one would also expect to be able to perform the 
equivalent calculations in the accelerated frame. It is this 
direction we will now pursue. 

We assume the transition to a higher energy state of 
the monopole, as in the photo-electric model, results in a 
current or cascade process that magnifies the occurrence 
of the quantum transition into a classical measurement. 
However, we don't model the quantum to classical pro- 
cess here. 

We assume the "binding energy" for this model is of the 
typical energy scale of the mode of the field we measure 
For the Gaussian in ( 10 1 we would have Ei — Eq 

The initial field state is taken to be the 
2| and the monopole is initialised into the ground state 
Pm{To) = \Eo){Eq\, so that the initial state of the sys- 
tem (monopole plus single field mode) is given by the 



§ 



Nhc 
L ■ 

mixed state 



product state p{tq) = Pm{To) (^ Tr_L^^|^)(\I'|. From the 
outset we work in the accelerated frame with the detec- 
tor positioned at ^ = ^q- In this frame the system will 
be assumed to undergo evolution according to the oper- 
ator of infinitesimal displacements of proper time. The 
different inertial and accelerating photo-count statistics 
will then be seen to arise because of the different proper 
time parameters under which the Hamiltonian systems 
evolve respectively. The probability of finding the sys- 
tem in some orthogonal state p/ = \E,x){E,x\ after an 
interaction time T can be calculated using a perturba- 
tion expansion of the equations of motion for the density 
matrix. The second order approximation is given by |_1.2]: 



V{^o,To;T) 






tq+T 



drji 



dri2y. 



(A3) 



{E, x|m(7yi)0(^o, Vi)p{Ta)m{V2)<i>{^o, V2)\E, x) + c.c, 

To this order the only processes that occur are those in 
which a single photon is either created or annihilated. In 
analogy with electrons excited in photo-diodes we assume 
a classical click of the detector occurs when the monopole 
makes a transition into any of the excited states. There- 
fore, to calculate the probability of a click we sum over 
all possible final energy states of the monopole. Further- 
more, since when we obtain a click we don't know which 
state the field has transitioned into the total single click 
probability is found by summing over all possible final 
states of the field. We then arrive at the probability for 
a single detector click in a time T given by: 



^(eo,ro;T) = |^^|(i?|m(0)|£;o)P-F(i?-i?o) 



(A4) 



where 



T{E - Eo) = 



To •'To 



(A5) 



Consider the terms in equation m\ since the accelerated 
worldline xq{t) = {S,o,t) lies completely in region I, the 
region II operators do not contribute to the monopole 



interactions (A2). Therefore, 

4'[x{t)] = / dkbkjWkji^o, t) + bljWkji^o, t)* = -0/(Co, t), 

(A6) 

and using the expansion: 



■>Pii^o,m)'4'i{^o,m) = :'0/(Co,?72)'/'7(Co,?7i) 



dk 

inlkl 



-i|fc|c(rj2-i)i) 



(A7) 



equation ( A5 1 can be rewritten: 

l-ra+T 



/■Ta + l pvi 

J To J Tn 



AE-Ea)(r,t-r)2)/h, 



To 



Tr 



{:^/(eo,^2)^/(eo,^i):Tri^ol*>(*l} 



dk 

47r|fc|' 



-i|fc|c(?)2-»)l 



+ C.C. 



(A8) 



It is useful to separate F{E — Eq) = Fn.o{E — Eq) + 
J^vac{E — Eq) into a normally ordered part and a vac- 
uum fluctuation contribution (that is independent of the 
state). First we focus on the vacuum fluctuation contri- 
bution. Deflning 77 = ?/i — rj2, we flnd: 



^vac{E — Eq 



dk 

47r|fc| 



I 

vacuum state, therefore: 

Tr{: ^Pi{^Q,V2)i'i{io,m) ■ Tr^.0ol*)(*l} 

+ ^d(Co, m)MCo,m)md^d\^) + C.C. (All) 
To analyse J-n.o{E—EQ) further, it is useful to decompose 



d?7(T- |77|)e 



\((E-Ea)/h+\k\c)ri t/ij-, into Rindler frequencies: 



(A9) 



Consider the innermost integrand. Since the time of 
measurement is much larger than the typical timescale 



of a transition, T ^ 



for all wave numbers k 



Ei—Eq ' 

the phase oscillates rapidly over the rj integration do- 
main. The inner integral is then approximately equal to 
2116 {{E — Eo)/H + \k\c). Since this is zero for all k the 
vacuum contribution is negligible. 

Now to evaluate the normally ordered field contribu- 
tion we use the expansion: 



0(e, r) = d>D(C, r) + Sr^i^, r) + ^' 



(AlO) 



In normal order all the annihilation operators are to the 
right of the creation operators, however, only the mode 
■01) is occupied in the state accessible to the detector, 
Trj^^j^|^)(vl'|. Therefore, no contributions to (A5l can 



arise from the 0' part of the field operator since those 
terms will be zeroed by acting on the orthogonal subspace 



V'd(C,^)= i^B, Wk,i)wk,ii^,T)dk. (A12) 



Consi der the contribution to (A8I from the first term in 
(|All|, let Auj = ^^ then: 



dm / rf^2e^(^-^°^^''^-''^^/'">D(^o, '72)^0(^0,^1) 

To -^To 

^ [ dkdk'{^B,Wk,l){^B,Wk'j) ^i(^^(k+k')-ro(\k\ + \k'\) 

4TTy/\k\\k'\ 






dr]2e 



-i(Auj+\k\)r]2 



(A13) 
Since the ip^ mode has a peaked frequency approxi- 
mately equal to the typical transitional frequency, uJc '^ 
Au, the innermost integral oscillates at double the op- 
tical frequency. But the measurement time is assumed 
to be many optical periods long (i.e., T ^ L/c = — ), 
therefore this term only makes a small contribution to 
the total integral. By analysing all the other terms in 
a similar way, it can be seen that the dominant contri- 



bution comes from the second term in (All I. One may 



recognise this step as the Rotating Wave Approximation, 
however this is a different approximation to the usual 
RWA in the inertial frame because the frequencies here 
are Rindler frequencies. One then obtains: 



J-(S - Eq) = (*|dtd|*) f dm r dme^''-''"^'^'''-'''^/^^U^Q,m + To)0D(Co,m + Tq) + C.C 

Jo Ja 



(A14) 



Note that all the dependence on the state of the field, |^), has completely factored out of the integral. The probability 
of a single click at the spacetime point (^0, '''o) (over a measurement duration T) then takes the form: 

ri^o,rQ;T) ^ a{^o,ro;T){^\d^d\^), 

where 



a(Co,T-o,r) 



M2 



Ek^I' 



\Eo) 



dm / dry2e^(^-^«)("^-''^'/V£(eo,r?2 + t-o)0d(6,'7i +To) 
Jo 



(A15) 
(A16) 
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We now consider two measurement scenarios. Consider 
first the case when the measurement time T is much 
shorter than the period of the entire pulse but still longer 
than the binding period, ^ ^^ <^ T <^ L/c. In this 
regime, the measurement time is too short to resolve the 
spectral line-width of the wave-packet and so the quasi- 
monochromatic approximation can be used: 



tpui^Q, m + To) ^ V'd(Co, T-o)e" 



(A17) 



where uJc = ^- In this case, a{^Q, tq, T) cx T. Therefore, 
the single particle absorption probability scales linearly 
with the measurement duration. 

The second case we consider is when the measurement 
time is taken to be greater than the time it takes the sin- 



gle mode to traverse the photo-detector, T > L/c. Then 
frequencies other than the peak frequency will contribute 



to the detection probability. In this case equation (Al5l 



must be interpreted as the probability of a single count 
per shot (unit time T). One can build up photo-count 
statistics by preparing an ensemble of identical systems 
and performing a single shot measurement (of time T) 
on each system in the ensemble. 

We conclude this appendix by noting that taking |$) 
to be the Minkowski vacuum state then agrees up to a 
device-dependent proportionality constant with the same 
expression arrived at using the generalised Glauber mode 
of detection ( 13 1 . 
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